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characterization of the invariant measures (in particular, all invariant phase separation 
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1. INTRODUCTION 

Cellular Automaton 184 (CA 184) is a discrete time process with state space 
{0, 1} Z and the following evolution rule: if r\ G {0, 1} Z is the state at time n then the state 
rj' at time n + 1 is defined by 

{1, if n(x) = r/(x + 1) = 1 
1, if rj(x) = 1 - r}{x - 1) = VxGZ (1.1) 

otherwise 

Here rj(x) denotes the value of?7:Z— >{0,l}at the coordinate x. The dynamics of CA 184 
will be denoted by the operator C : {0, 1} Z — > {0, 1} Z defined by Crj = rj' . Thus, formally, 
CA 184 refers to a sequence {q n , n G N} such that n n+ i = Crj n , Vn. 

CA 184 models deterministic motions of identical particles on Z that obey the following 
rules: there may be at most one particle per site, and at each integer time each particle 
inquires whether the site of Z to the right of its current position is empty of another 
particle, and if it is so then it instantaneously jumps to this site. These rules follow 
immediately from (1.1), if one interprets r/(x) = 1/0 as " presence/absence of a particle at 
the site x G Z in rf J . 

The number "184" in the name of this process is due to the classification of Wolfram 
[Wl] (see also [W2]) of a class of cellular automata. CA 184 has been used to model traffic 
([N], [Nl] and references therein) and to classifies densities in binary strings ([CST]). CA 
184 has two "stochastic" counterparts. One is the Totally Asymmetric Simple Exclusion 
Process (TASEP). The second one has the same evolution rule as CA 184 but with "a noise" 
that is introduced by setting that each particle that can jump will do so with probability 



f Instituto de Matematica e Estatfstica, Universidade de Sao Paulo, Cidade Universi- 
taria, Rua de Mataao, 1010, 05508-900, SP, Brazil, pablo@ime.usp.br, vbell@ime.usp.br 



1 




p independently of anything else; we thus call this process "CA 184 with noise" (CA&N). 
The invariant measures for TASEP and CA&N have been characterized in [L] and [Y], 
respectively (curiously, only the case p < 1/2 was studied in [Y], but we believe that 
a similar technique allows to extend the results to any p) but the approach employed 
there (which is a stochastic coupling) does not apply to CA 184 (because of the lack of 
stochasticity in its dynamics). 

Annihilating particle system studied here can be also found in the literature under 
the names Ballistic Annihilation (BA) and Annihilating Deterministic Motions. We shall 
adopt here the name BA. BA is a discrete time process with state space {—1, 0, 1} Z and 
the following dynamics: if C is the state at time n then the state (' at time n + 1 satisfies 

if C(x — 1) = 1 and neither ((x) = —1 nor both ((x) — 0, ((x + 1) = — 1 
if ((x + 1) = — 1 and neither ((x) = 1 nor both ((x) — 0, ((x — 1) = 1 
otherwise 

(1-2) 

for all x G Z. Here ((x) denotes the value of ( : Z — > {—1, 0, 1} at the coordinate x. The 
dynamics of BA will be denoted by the operator A : { — 1, 0, 1} Z — > { — 1, 0, 1} Z defined by 
AQ = Thus, formally, BA refers to a sequence {( n , n G N} such that Cn+i = A( n , Wi. 

BA may be also interpreted in terms of particles. We shall call them A-particles in 
order to distinguish them from those that move in CA 184. The values 0, 1, —1 of ((x) 
are interpreted by saying that the site x is respectively, free of an A-particle, contains an 
A-particle with velocity 1, and contains an A-particle with velocity —1. In the terms of 
particles, the dynamics of BA acquires the following interpretation: each A-particle moves 
along R with its velocity (going in the direction to — oo (+oo), if the velocity is negative 
(positive, resp.)) and annihilates when meets another A-particle; upon annihilation both 
A-particles disappear from the system forever. To link this interpretation with (1.2), it is 
necessary to note that if A-particles occupied only the sites of Z then after a unit of time 
has elapsed the survived A-particles are again solely at the sites of Z. 

BA is a natural model for the reaction A + B — > inert. For an extensive list of 
references in this respect, we refer a reader to [EF]. A modification of BA in which the 
set of possible particle velocities is larger than {—1,1} has been considered in the works 
[BRL] , [KSL] . We also mention here the work [E] which studies a coalescing particle system 
employing a very simple relation between it and BA. This coalescing particle system is a 
process in which particles move as in BA but coalesce instead of annihilate, upon collision; 
coalesced particles choose then a new velocity from {—1, 1} with equal probabilities. 

Surface Growth Model (SG) is a discrete time process with state space 1Z, the 
space of piecewise linear functions (from R to R) that have slope either +1, —1, or 
between any two consequent integer abscissas, and the following dynamics: If /(•) is the 
time n then /(•) G 71, the state at time n + 1, is obtained from /(•) by the following 
rule: if x G Z is a local minimum of /, then the values of / at [x — 1, x + 1] are obtained by 
reflecting the graph of / around the straight line which links the points (x — 1, f(x — 1)) 
and (x + 1, f(x + 1)); at any point y G R which does not have a local minimum of / 
in its neighborhood of radius 1, we set f(y) = f(y). SG will be the name of a sequence 
{f n (-), n G N} such that / n +i(-) = /„(•) , Vn. 
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Imagine /(•) as the surface of a two-dimensional solid above some reference horizontal 
line. Imagine then that diamond shaped particles of the side length 1 are thrown on this 
solid, and those of them that fall in local minima of /(•) stick to the solid while others 
disappear. The surface of the new solid will be then what we have defined as /(•). This 
justifies the name of the process just introduced. We observe that it is also known under 
the name polynuclear growth model (PNG) and (one of its modifications) was first studied 
almost 30 years ago in [Ber] (see also [KS2] for a list of more recent works). 

Equivalence of dynamics of two processes means here an existence of a transforma- 
tion T from the state space of one process to that of another one such that if Tr\ = £ then 
Tvj = where rj' and £' are the states of these processes at time n + 1 given their states 
at time n were rj and ( respectively. The fact that the dynamics of CA 184, BA and SG 
are equivalent has been known since at least the work [KS]. We present it in Assertions 1.1 
and 1.2 below. The proofs are straightforward and thus, omitted. 

Assertion 1.1. (equivalence between CA 184 and BA) Define Ti 84 ba '■ {0, 1} Z — > 
{-1,0,1} Z by 

{Tl8A,BAV)(i) = l-7/(i)-7/(i+l), zGZ (1.3) 

Then, if for some rj and (, Ti 84:j ba(v) = C then Ti 8 4^ba(v') = C> where rj' relates to rj via 
(1.1), and (' relates to ( via (1.2). 

Remark 1.1. Tis4,ba is an injection but not a surjection. Indeed, if any two conse- 
quent sites of Z are occupied by A-particles with opposite velocities in a configuration 
( G {—1,0, 1} Z then there does not exist a rj G {0, 1} Z such that Ti 8 ^baV = C- As a 
consequence, a certain care is necessary, when CA 184 is studied with the help of BA (the 
reader will see it in the proofs of Theorems 2.2 and 2.3). 

Assertion 1.2. (equivalence between BA and SG) Define T B a,sg '■ { — 1, 0, 1} Z — > 1Z of 
the following form 

(T B A,SG()(n) - (T BA ,SG0(n - 1) = CH, n G Z 
(T B a,sg()(0)=0 

(certainly, it is sufficient to determine f G 1Z on Z). Then, if Tba,sgC = / + a f° r some 
( G {-1, 0, 1} Z , / G K and a G R, then T BA ,sgC = f + P, where (' relates to ( via (1.2), 
f relates to f via the dynamics of SG and f3 is a real number whose value depends on 
a and f .(Above, the relation g = h + a for functions g and h and a number a means 
g(x) = h(x) + a\fx G RJ 

Remark 1.2. The actual values of the constants a and (3 are ignored because they will 
be of no importance, when SG is employed to investigate BA and CA 184, and also when 
BA or CA 184 is employed to investigate phenomena related to modification of the shape 
of the surface in SG. In this respect we observe that in SG, the surface both grows and 
changes its shape, but only the phenomena related to the shape change will be of interest 
to us. 

The contents of this paper: Section 2 characterizes the invariant measures for the 
dynamics of the processes defined above. This is done with a help of a quite elementary 
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property of BA presented in Assertion 2.1; it states at any time n the distance between 
two diverging subsequent particles (i.e., no any other particle in between) cannot be less 
than 2n + 1. (We note that this property has been noticed in [CST] where it was employed 
to show that CA 184 on a ring of size N will come to its invariant state by time not greater 
than N/2.) The set of the invariant measures which are also translation invariant is then 
easily classified (Remark 2.2). For BA, an extremal measure of this set is supported either 
by configurations that have only particles with velocity +1 or by those that have only 
particles with velocity —1, moreover, under such a measure the positions of these particles 
form a stationary ergodic process on Z. Applying then Assertions 1.1 and 1.2 to this result 
one easily characterizes the translation invariant measures that are invariant for CA 184 
and SG. For extremal translation invariant measures of CA 184 we calculate the flux of 
particles through zero. The curve of the flux (as a function of the density of particles 
p G [0, 1]) is 1/2 - |l/2 - p\; we observe that it is known to be p(l - p) for TASEP. 

The invariant measures of BA which are not translation invariant are called "phase sep- 
aration" measures. These measures are supported by configurations in which all particles 
with velocity +1 are situated to the left of all those whose velocity is —1. Our characteriza- 
tion of the phase separation measures is obtained via introducing in BA an extra "second 
class" particle in the way such that its paths record the positions of particles in BA, and 
studying the properties of these paths. The second class particle is always between the 
rightmost particle with velocity +1 and the leftmost particle with velocity —1. It moves to- 
gether with the latter till it is annihilated by the former. At this moment, the second class 
particle changes its velocity and starts moving to the right until it meets another particle 
with velocity — 1. The definition of the second class particle yields the following relations 
of its path to particle positions in BA: each time interval the second class particle moves to 
the right (respectively, left), is equal to the distance between two subsequent particles with 
velocity —1 (respectively, +1). As for the relative positions of particles with velocity —1 to 
those with velocity +1, it must be such that the their annihilation times form a point pro- 
cess identical to that formed by the times of change of velocity of the second class particle. 
These relations are employed in Theorem 2.3 which states that BA is distributed due to a 
phase separation invariant measure if and only if the second class particle paths belong to 
the set E := {e = (ei)i e z : e 2 % G 2Z Vz G Z, and — ej+i G {—1, 1} Vz G Z} and have the 
distribution which is invariant with respect to the shift of E by 2 (i.e. — > e^i). Above 
ej is interpreted as the position of the second class particle at time i. For CA 184, a similar 
results holds with the only difference that now the invariance is with respect to shift by 1, 
and E := {e = (e^igz : e.% — G {—1, 1} Vz G Z, and if for some i and j, j > i it holds 
that e; - ei-i = -(e i+1 - e*) = e i+2 - e i+1 = . . . = ej - e_,-_i = -(e J+ i - e,) then j - % 
is odd}, other words, the second class particle motion in CA 184 is a time-homogeneous 
process with the state space E. We note that the ideas similar to those just presented, 
have been used in [FR] to construct phase separation measures for a so-called Boghosian 
Levermore cellular automaton. 

In Section 3 we review and a generalize results from [BF]. Theorem 3.1 presents a 
hydrodynamic limit for the shape of surface in SG. It says that if in SG {f n {')i n G N} 
the distribution of c n /o(n-) converges to some process W(-) for some sequence of numbers 
{c n } then c n f*(n-) converges to the process M /mm (-) which is defined by W rnin (x) := 
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min{W(y) , x—1 < y < x+1}, Vx G R; above /* denotes a particular function whose shape 
coincides with the shape of f n . The cornerstone of this result is the following property 
of SG: if for a function /(•) G 1Z we define the function g(-) by g(x) = min{/(y) , x — 
1 < y < x + 1} then g(-) and /(•) have the same shape. With help of Assertions 1.1, 
1.2 we then derive from Theorem 3.1 the time — > oo limits for the particle distribution 
in BA (Theorem 3.2) and in CA 184 (Theorem 3.3), under appropriate rescaling (i.e. 
hydrodynamic limits). These theorems express the limit of distributions of particles in 
BA and CA 184 in terms of a limit of their counting processes. A counting process for 
BA {Cn,n G N} is, by definition, the process {Tba,sg(Cti), n G N}. Since the latter 
closely relates to SG, as specified in Assertion 1.2 and Remark 1.2, then Theorem 3.2 is 
immediate from Theorem 3.1. To obtain a counterpart for CA 184, just a simple trick is 
required. The results of Section 3 have two natural applications which we demonstrate via 
examples. Example 3.1 shows how BA and CA 184 may be employed in order to find the 
law of VF mm (-) for a given process W(-). Example 3.2 demonstrates a simple way to find 
approximately the particle distribution in BA or CA 184 for large times. 

In Section 4 we calculate the rate of convergence of CA 184 starting from the Bernoulli 
1/2 product measure, to its invariant state. This is equivalent to calculating the rate of 
the decay of particles in BA with a particular initial state. The argument we employ is not 
new. It has been used to estimate analogous rates in various models (see [Fi], [KS], [EF]). 
This argument is based on estimating a certain quantity which we specify in Remark 4.2. 
We observe (in Remark 4.1) that the result discussed in this section could be also obtained 
using the tools presented in Section 3. 

We consider the results of Theorems 2.3, 2.4, 3.3 to be our novel contribution to the 
study of the dynamics of BA, CA 184 and SG. Theorems 2.1, 2.2 and 3.1, 3.2 have original 
form but the bulk of their proofs is based on ideas and techniques that have appeared in 
literature related to the field. 

2. INVARIANT MEASURES 

In this section, we give the necessary and sufficient condition for a measure to be 
invariant for BA (Theorem 2.1) and for CA 184 (Theorem 2.2); the corresponding result for 
the shape in SG may be easily obtained from that for BA (by using T from Assertion 1.2) 
and thus, will not be presented here. Those invariant measures which are translation 
invariant are characterized in Remark 2.2, and those which are not translation invariant 
are characterized in Theorems 2.3 and 2.4. At the end of Section 2 we calculate the flux 
of particles in CA 184 when it is distributed with an invariant and translation invariant 
measure. 

We start with some definitions. Elements from {0, 1} Z and from {—1, 0, 1} Z will be 
called configurations, and their values at sites of Z will be interpreted in terms of particles as 
indicated in Introduction. A-particles of BA will be called simply particles, and those whose 
velocity is +1 (resp., —1) will be called positive (resp., negative). We set © := {—1, 0, 1} Z . 
By X we will denote the set of the measures that are invariant for the dynamics of BA. 
To say fx G X means that / u(A~ 1 (L r )) = fi(ll) for any cylinder set U C (that is, a set 
of the form U := {C G 6 : C(*i) = a i> • • • > C(*fc) = for some k G N, i\, 12, . . . , %k G Z 
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and ai,a2, ■•• ,a,k G {—1,0,1}.) By t\ we will denote the translation of Z by 1 to the 
right. By X we shall denote the set of those measures on {-1,0, 1} Z that are translation 
invariant, i.e., h(ti 1 U) = fi(U) for any cylinder set U C 0. By we shall denote the 
configuration that does not have any particle: = OVz G Z. By + := {0, 1} Z \ {</>} 

(resp., 6_ := {— 1,0} Z \ {</>}) we denote the set of those configurations of which have 
solely positive (respectively, negative) particles. By V+ (resp., V-) we denote the set of 
measures concentrated on 0_|_ (resp., on 0_). Given a configuration (6 0, we will call 
two particles of £ subsequent, if £ does not have any other particle between these two. A 
pair of subsequent particles is said to be converging if the leftmost particle of the pair is 
positive and the rightmost one is negative. The diverging type is introduced in analogous 
way. By S we denote the set of those configurations (66 which satisfy the following 
two conditions 

(i) there is only one pair of subsequent converging particles; 

(ii) both the number of positive particles and the number of negative ones in ( are infinite. 

By S we will denote the set of measures concentrated on S . These measures are custom- 
arily called phase separation measures. 

Theorem 2.1. fx G X if and only if for some 011,0:2, 0:3, «4 > such that 01+02+03+04 = 
1 it holds that 

^ = 01^1 + 02^2+03^3 + 04^, (2.1) 
whereas fi\ G V+ fl X, ^ 2 eP_nT, ^3 G S n X n T c . 

Remark 2.1. T c means the complement to T. We wrote S C\lr\T c in order to emphasize 
that 5fll cannot contain a translation invariant measure. This fact can be established 
by the following reasoning: Let pos(Q and neg(() denote the positions of respectively, the 
positive and the negative particles of the unique pair of converging particles from (6 s . 
Define the function I:O s ^ZU(Z + 1/2) by X(£) := (neg(C) - pos{Q)/2. Now, if ( 
is distributed by a translation invariant measure then the distribution of X(Q would be 
invariant with respect to translations of Z, which is impossible. 

Remark 2.2. The set (V+ fl T) U 5^ (and (V- fl X) U 5^, by analogy) is by its definition, 
the set of all translation invariant measures on {0, 1} Z . The extremal points of this set is 
known to coincide with the set of all ergodic (with respect to T\) measures on {0, 1} Z . 

Proof of Theorem 2.1. Observe that 

if CG0+thenri(C)=A(C) (2.2) 

It follows then easily from (2.2) that if \i\ G P+flX then \i\ G X. Analogously, if \i2 G "P-flX 
then fi2 G X. Since ^3 in (2.1) is picked from a subset of X and since obviously 5^ G X 
then fj, defined by (2.1) must belong to X. This completes the proof of the "if" part of 
the theorem. The "only if" part is more difficult. It will be proved with the aid of the 
following 

Assertion 2.1. Let Co G and m G N be arbitrarily fixed. Let ( m denote the config- 
uration of particles in the BA at time m, starting from ( , i.e., Cm = ^. m Co- Then the 
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distance between the particles of any pair of diverging particles from ( m is not less than 
2m + 1. 

Proof of Assertion 2.1. Find any pair of diverging particles in ( m . Observe that the initial 
position of the negative particle of this pair must be to the left of that of the positive one, 
because if not then these particle would have annihilated each other by time m. The 
assertion follows from this observation and from the fact that these particles diverge with 
velocity 2. X 

Continuation of the proof of Theorem 2.1. Assertion 2.1 provides that if /i G I then it 
gives weight zero to any configuration that has at least one pair of diverging particles. 
Thus such ii must be concentrated on © + U 0_ U6 S U {</>}. Since the dynamics A does 
not mix these sets (that is, there is no a configuration £ from one of these sets such that 
A( belongs to another one) then the expansion (2.1) follows. X 

Let T be defined by (1.3). We define A C © as the image of {0, 1} Z by T. Namely, 
A := {( G 6 : 3r/ G {0, 1} Z such that Trj = £}. ^From this definition and that of T, A is 
easily characterized: A := {( G 6 : the distance between any subsequent particles with the 
same velocity is odd, and the distance between any two subsequent particles with opposite 
velocities is even}. The set A is used in the following statement: 

Theorem 2.2. A measure \i on {0, 1} Z is invariant for CA 184 if and only if 

H = ai/ii +ct2H2 + «3/«3 + a4 ( - 77 e ) , for some cti, «2, «3, 014 > 0, ot\ +ct2 + «3 + a4 = 1 



where the measures [i\ and are translation invariant and are concentrated on the sets 
T _1 (0 + fl A) and T _1 (6_ fl A) respectively; the measure ^3 is invariant for CA 184 
but not translation invariant and is concentrated on T _1 (G S fl A); the measure 5 Q (resp., 
5 e ) gives mass 1 to the configuration o (resp., e) where o and e G {0, 1} Z are such that 
0(1) = 1, e(i) = 0, when i is odd, and o(i) = 0, e(i) = 1 when i is even, i G Z. 

The proof of Theorem 2.2 comes out by a reformulation of the proof of Theorem 2.1 in 
terms of CA 184. This reformulation is based on the property of T stated in Assertion 1.1, 
and the facts that T({0, 1} Z ) = A and that T is invertible on {0, 1} Z \ {o,e}. The only 
care thus, should be taken because of the fact that T(o) = T(e) = 0. But this causes no 
difficulties in the mentioned reformulation because we have that C(o) = e and C(e) = o 
and as a consequence of this, we also have that if n is invariant for CA 184 then /x(o) = //(e). 
* 

We now characterize the set of invariant measures of BA which are not translation 
invariant. Set 

E:={e = (e^ezuz+i : e, G Z Vz G Z and e % - e l+ i G ^} Vz G Z U Z + ^} 

and define the shift te '■ E — > E by TE(e) = e' e\ = e^+i Vi. Denote then by £ T the set 
of the measures on E which are invariant with respect to the shift te- 
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Theorem 2.3. There is a bijection F* : S PI J — > £ T . The mapping F* will be constructed 
explicitly in the proof of the theorem. 

Proof. The idea is to introduce a second class particle in BA and to show that BA 
is distributed by an invariant phase separation measure if and only if the second class 
particle's motion has invariant distribution with respect to shifts of time. The second class 
particle in BA is defined by the following rules: (%) it is an extra particle that does not 
affect the evolution of other particles, and (ii) it moves with the velocity that changes 
from +1 to —1 and back depending on the BA particles that the second class particle 
meets, namely: (ii-a) when the second class particle meets a negative particle it changes 
its velocity to — 1 and keeps going together with the met particle; (ii-b) at the time this 
met particle is annihilated, the second class particle changes its velocity to +1. One more 
postulate is: (in) if two particles annihilate at a time n G Z U Z + \ then at this time, the 
second class particle is at the annihilation point and changes its velocity from —1 to +1. 
We want the rules (i)-(iii) to define correctly and uniquely the second class particle and 
moreover, we want each trajectory of the second class particle to determine uniquely the 
configuration of BA at any time. To achieve these objectives we have to consider BA with 
double-infinite time and to record the second class particle position at times from the set 
Z U Z + |. Note that Z + | has to be added because annihilations may occur at times 
from this set. Also note that the distribution of the second class particle trajectories turns 
out to be invariant with respect to the shift of time by 1 but not necessarily to the shift 
of time by \. The details are presented below. 

We start introducing a terminology which we shall need to construct BA with double 
infinite time. For an element u from {— 1,0, l} ZxZ , let us write u n {m) for its values at 
(n, m) G Z x Z. Let us also write u n for the restriction of u G {— 1,0, l} ZxZ obtained 
by fixing its first coordinate to n; observe u n G { — 1,0, 1} Z . Define then ft := {u G 
{— 1, 0, l} ZxZ : uj n G S and w n +i = Au n Vn}. Fix now an arbitrary [i G S flX. There is 
a standard procedure to construct a probability measure (on an appropriate a-algebra 
T of subsets of ft) with the following properties (2.3), (2.4): 

P£ [u : u n G D] = [i [D] , Vn G Z and V cylinder D C {-1, 0, 1} Z (2.3) 

If we define the shift operator tq, : ft — > ft acting by (rn(o;)) n (m) := u n -i(m), Vn, m, then 

P£[£>] =P£[rn 1 (£>)], V cylinder D C { — 1, 0, l} ZxZ (2.4) 

P(^ will be called BA with double infinite time and the marginal distributions fx. Such a 
choice of the name becomes obvious if one interprets a; as a trajectory of a certain process, 
and u> n (m) as its value at time n and the location m. We recall that it is known that given 
[i G iSflZ, there is a unique P^ satisfying (2.3), (2.4), and, moreover, given any probability 
measure P on ft which is invariant with respect to rn, there is a unique \i G S D X such 
that P(^ = P, other words, n — > P^ is a bijection from S fl X to the set of all probability 
measures on ft which are invariant in respect to tq. 

By certain technical reasons that will become clear soon, we need to have a record of 
the positions of particles in BA also at times from the set Z + |. We thus associate to 



8 



each wGfian element uj' E { — 1, 0, l} ZuZ +2 xZuZ +2 by the following rule: for each n E Z 
we set 

, / w n (m), when m E Z 

^(™) = U whenmGZ + I ^ 



and 



if m E Z 



/ , s _ I 1, when mGZ + i and if cu n (m — |) = 1 and uj n (m + \) ^ — 1 , 
-1, when m E Z + \ and if w n (m + |) = -1 and u n (m - \) ^ 1 ^ J 
1 0, otherwise 

The relation between uj and uj' may be explained in the following manner: Suppose that two 
persons observe the evolution of the same process BA. Suppose that one of them records 
the states of the process at times Z while the other one records the states of the process 
at times Z U Z + |. Then, uj and uj' are the histories recorded by respectively, the first and 
the second person. Having this interpretation in mind and observing that a configuration 
of particles in BA at time n E Z U Z + | determines uniquely the configuration at time 
n + \ we conclude that the mapping uj — > uj' is a bijection. 

We are now in a position to define the mapping F : O — > F. To obtain e = F(uj), we 
construct first a;' by the rules (2.5)-(2.6) and then, for each neZUZ + |,we define e n 
by the following rules: 

1 if a point rnGZUZ + ^isan annihilation place of two particles at time n in uj' (that 
is, if uj' n (m) = 0, uj' n _ 1 (m — |) = 1, ^j' n _\{jn + \) = —1) then we set e n = m; 

2 if e„ i = m and u/ x (m) = —1 then we set e n = e„ i — 4; 



3 if e„ i = m and a;' i (m) = then we set e n = e„ i + i. 



An important observation is that F is a bijection. This follows from two facts. The first one 
is that uj — > a;' is a bijection, as mentioned above. The second fact is that uj' determines 
uniquely e E E and vice versa. This fact may be verified straightforward. We observe 
that in order that F~ 1 e be uniquely determined, it was essential that the elements of E be 
indexed by the set ZUZ+^. (We suggest a reader to choose an arbitrary e E E, to construct 
then e E E by changing the value of e at |, and to verify then that F~ 1 (e) ^ F~ 1 (e).) The 
mapping F introduces naturally a measure on E via [£)] := P^ , V cylinder 

D C E. But because of the property 

^(F-^e)) =F- 1 (r £ (e)) 

that follows from our constructions and because of (2.4), we conclude that P^ [£)] = 
V^ E [ T - 1 D]^D E E, that is P£ is from E T . We thus define F* by F*(ji) := P^, for 
(j, E S fl X. The fact that F* is a bijection follows from that the mappings F and /i — > P^ 
are bijections. X 

Recall Theorem 2.2. Let us give the name invariant phase separation measure of 
CA 184 to an invariant measure of CA 184 which is concentrated on T _1 (0 S fl A). The 
following result characterizes these measures. 
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Theorem 2.4. There is a one-to-one mapping between the set of all invariant phase 
separation measures of CA 184 and the set of all stationary processes (X n ) ne z with values 
in Z which satisfy the conditions 

(a) X n+ i - X n e {-1, 1}, for all n; 

(b) if for some j > i, it holds that Xi — Xi-\ = — (X i+ i — Xi) = X i+2 — X i+ i = . . . = 
Xj — Xj-i = —(Xj+i — Xj) then j — i is odd. 

This mapping can be obtained as an appropriate restriction of F* in the manner outlined 
in the proof. 

Proof. The mapping T from (1.3) induces a natural mapping between the distributions of 
CA 184 and that of BA. The property that T : {0, 1} Z — > A is a bijection (recall, A has 
been defined before the statement of Theorem 2.2) then yields that the induced mapping 
is a bijection between the set of all invariant phase separation measures of CA 184 and a 
subset of (S fl J) which we shall denote by (S fl X) 184 . 

Thus, Theorem 2.4 is proved, whence we show that F*((S fl I)i8a) is the set of the 
processes X as defined in the theorem's statement. This can be done directly via repeating 
the proof of Theorem 2.3, but applying it to (S flX)i84 instead of S flX. For this application 
to be correct, it is necessary to note the difference between these two sets. It follows from 
the property of T mentioned in the above paragraph that a measure from S fl X is in 
(S H X)i84 if and only if it gives weight to the set A c = {( G {— 1, 0, 1} Z : the distance 
between any subsequent particles with the same velocity is even or the distance between 
any two subsequent particles with opposite velocities is odd}. This suggests the following 
two modifications in construction of F*((«S nX)i 8 4) compared to that of F*(S(~)l): First, 
we do not need to record the second class particle positions at times Z + \, because there 
are no annihilations at these times. Second, the set of paths of the second class particle 
should be modified so that when a path from this set is translated into corresponding 
configuration of annihilation particles, then the latter must belong to A c . To satisfy these 
requirements, we imposed the conditions (a) and (b) on the process X. The rest of the 
proof goes almost identically to that of Theorem 2.3 and thus, will be omitted. A 

Remark 2.3. We recall (from [L] and [Y], respectively) that in the case of TASEP and 
CA&N, the set of phase separation measures is {i>^ n \ — oo < n < oo}, where gives 
mass 1 to the configuration such that rj^ n \x) = lVi > n and rj^(x) = OVx < n. 

Flux of particles in CA 184 which is distributed due to an translation invariant measure 
fi which is invariant for CA 184, is defined as 

^(lXl-Tfc-iCl))} 

We shall now show that the flux of particles in CA 184 is given by 1/2 — 1 1/2 — p\, if 
CA 184 distributed by such a measure p, where p denotes the density of particles. (Here 
the density of particles is understood as p{i] n (i)} which certainly does not depend on n and 
i.) Let p, be such a measure and consider first the case p < 1/2. Then Theorem 2.2 implies 
that p £ T _1 (0_ fl A) which in turn, implies that p{n n (i) = | i] n (i — 1) = 1} = 1, Vz Vn. 
Thus the flux of particles is p. In the case when p > 1/2, the flux of particles is 1 — p. This 
results from the following reasoning: We define the flux of holes as p{(l — i] n (0))ri n -i(0)} 
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and derive that it equals 1 — p, this result follows by adapting the reasoning that have been 
used for the case p < 1/2. We then observe that the flux of holes must be equal to the flux 
of particles (this observation may be justified in one of the following ways: either using 
the fact that p is invariant, or from a direct comparison of the definitions of the fluxes). 
It is left to consider the case p = 1/2 which occurs as follows from Theorem 2.2, when 
p = l/2(5 + S e ). When it happens, then flux is obviously 1/2. 

3. HYDRODYNAMIC LIMIT. 

This section concerns the limit (as time — > oo) shape of surface in SG (Theorem 3.1) 
and the limit distribution of particles in BA (Theorem 3.2) and in CA 184 (Theorem 3.3). 
These results are called hydrodynamic limits because along with time tending to oo we 
apply a certain rescaling that depends on time, to the quantities whose temporal limit is 
of interest. We also demonstrate here applications of these results (Examples 3.1, 3.2). 

For y > 0, let us denote by M y the operator that brings a function /(•) : R — > R to 
the function g(-) : R — > R in the way such that 

g(x) = (M y f) (x) = min{ / (z) : x - y < z < x + y}, VxeR 

Let us then introduce Surface Modification Process (SM) as the name of a sequence 
{fn{')i n £ N} of random functions from 1Z such that / n +i(-) = (Mi/ n )(-) for each n G N. 
One then can check straightforwardly the following property: 

for any f e 71 there is an a G R such that (Mi/)(-) = /(•) + a (3.1) 

whereas the construction rule of /(•) has been specified in the part of Introduction that 
defined the Surface growth process. The above (3.1) says that at each time, the shape of 
the states of SG and SM is the same, provided they stated from the same state. However, 
SM is more handy for formulating the hydrodynamic limit for this shape. To state this 
result, we need first to introduce some notations: (1) For a process W(x), x G R, we 
will denote by W min (-) the process defined by W min (x) := (M 1 W)(x), \/x G R. (2) The 
convergence of processes discussed in this section and throughout the paper, is understood 
as the weak convergence of their distributions on C[a, b], the set of continuous functions on 
[a, b], for any — oo < a < b < oo (see [Bi]). 

Theorem 3.1 (hydrodynamic limit for the shape of surface in SG). Let /o(-) be a random 
function with the state space 1Z such that c n f (n-) — > W(-), as n — > oo, for some sequence 
{c n , n G N} of real numbers and some process W(x), i6R. Then the SM {f n (-), n G N} 
satisfies the following property 

c n fn(n-) - V^ min (.) (3.2) 



Proof. First, we note that for any function h : R — > R, the modulus of continuity (which 
definition one finds in [Bi]) of M\h does not exceed that of h; this property of Mi can 
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be verified straightforward. This fact allows us to apply Theorem 5.1 from [Bi] to the 
assumption c n f (n-) — > W(-) to conclude that 



We must comment on the usage of two " •" in (3.3) in order to avoid a possible confusion 
in its interpretation: Mi(c„/o(n-))(-) means the function obtained by rescaling f by c n 
along the axis of ordinates and by n _1 along the axis of abscissas, and then, by applying 
Mi. It may be verified directly that for any continuous function h and any constant c, 



and also that n subsequent applications of Mi is equivalent to one application of M n . 
Thus, (3.2) follows from (3.3). * 

Let (O, JF, P) be an abstract probability space. We shall say that a random function 
f n (-) : O — > 1Z is a counting profile of Cn : ^ — ► { — 1,0, 1} Z (resp., rj n : O — > {0, 1} Z ) if 



/(*)[<*;] -/(fc-l)[u;] = C(fc)M (resp., /(fc)[u;] -/(fc-l)[u;] = -rj(k)[<jj]), V/c G Z (certainly, 



it is enough to specify the values of a function from TZ just on Z) and Vcu G O, where 
we denoted by Y[oj] the realization on uj G O of a function y defined on We say that 
{fn(-), n G N} is a counting process of BA {C n , ft G N} (resp., CA 184 {i] n , n G N}) if for 
each n, f n {-) is a counting profile of ( n (resp., rj n ). 

^From the definitions just introduced, Assertion 1.2 and (3.1), we conclude that if a 
random function /o(-) is a counting profile of a random configuration Co then SM {/n(-), n G 
N} is a counting process for BA {Cn, G N}. Thus, Theorem 3.1 leads to the following 

Theorem 3.2 (hydrodynamic limit for particle distribution in BA). Let {Cn, n G N} be a 
BA. Assume there is fo(-), a counting profile of Co, such that c n fo(n-) — > W(-) as n — > oo, 
for some sequence {c n , n G N} of real numbers and some process W(x),x G R. Then 
there is a counting process {/ n (-), n e ^} of BA such that c n f n (n-) — > VF mm (-). Precisely 
to say, this counting process is SM with the initial state fo(-). 

The next result is a counterpart of Theorem 3.2 for CA 184. Note that it contains 
an additional assumption "c n — > as n — > oo". In respect to this note, we observe that 
in all applications of hydrodynamic type limit theorems we are aware of, this condition is 
satisfied. 

Theorem 3.3 (hydrodynamic limit for particle distribution in CA 184). Let {rj n , n G N} 
be a CA 184. Assume there is go(-), a counting profile ofrjo, such that c n go(n-) — > W(-) as 
n — > oo, for some process W(x), x G R, and some sequence {c n , n G N} of real numbers 
satisfying c n — > as n — > oo. Then there is a counting process {<7n(-)> n e ^} °^ ^4 
such that c n g n {n-) — > VF mm (-). The construction rule of g n (-) will be specified in the proof 
of this theorem. 

Proof. It is straightforward to see that 



M 1 (c n f (n- 



))(.)^( Ml W)(-) = W mi \) 



(3.3) 



(Mic/i(n-)) (x) = c(M n h) (nx), Vx G R, 



if T?] = C for some n G {0, 1} Z and C G {-1, 0, 1} Z , where T is from (1.3), 
and if g and / are counting profiles of rj and C respectively, 
satisfying g(0) = f(0) then \ f(x) - g(x)\ < 1 Vx G R 



(3.4) 
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Let (o := Ti] , where T is from (1.3), and let /o(-) be the counting profile of Co such 
that /o(0)[u;] = go(0)[u>],Vu>. By theorem's assumptions and (3.4), c n fo(n-) — > W(-). The 
latter and Theorem 3.2 yield that c n f n (n-) — > W min (-) for an appropriate counting process 
{fn(-), n £ N} of the BA {C n , n e N}. This conclusion, (3.4) and the assumption c n — > 
then yield that c n g n (n-) — > VF mm (-), if we define <7 n (-) to be the counting profile of r\ n such 
that <7n(0)M = /„(0)H Vw. * 

Two applications of the above results are quite natural. We shall demonstrate them 
in examples below. 

Example 3.1 demonstrates that given a process W(-), the distribution of VF mm (-) can be 
found by considering an appropriate cellular automaton. This example is a brief review 
of a result contained in [BF] in which W(-) was taken to be B(-), the one-dimensional 
Brownian motion with B(Q) = 0. 

Let Cn-,n G N, be the BA whose initial state has the following distribution: 

Co(i),ieZ, arei.i.d. withP[CoW = 1] =P[CoW = -1] = V 2 ( 3 - 6 ) 

It is then not hard to find exactly the distribution of Q n = A n Q G . It is as follows: it is 
invariant with respect to translations and reflections of Z, and given a site % is occupied 
by a particle, the probability that the nearest particle to its right has the same (opposite) 
velocity is (1 + U2 n )~ 1 (resp., U2 n (l + U2 n )~ 1 ) and can be expressed exactly as a function 
of U2n, where v>2 n is the probability that a simple symmetric one dimensional random 
walk returns to the origin at time n; this expression will depend on the velocities of the 
particles under the consideration. For details we refer a reader to Theorem 1 in [BF] . Two 
fact contributed to the success in finding the law of ( n in such an exact form. First, a 
simple expression for the distance between a given particle and its annihilating companion 
(this is the essence of Assertion 4.2 presented in Section 4). Second, the relation of the 
distribution of this distance to ii2 n which stems from a particular choice of the law of Co, 
(3.6). 

Let now / n (-) denote the counting profile of Q n that satisfies / n (0) = 0. Clearly, 
the law of / n ( - ) can be derived from that of ( n presented above. It turns out that the 
known asymptotics for U2 n (see [F]) is sufficiently precise to enable one to derive then the 
distribution of the process to which n _1 / 2 / n (n-) converges (for details, see Theorem 2 in 
[BF]). If f n were M n f Q then this law would be of £? mm (-), by Theorem 3.2 and because 
n -1 / 2 /o(^') —> B(-). By technical reasons however (recall Remark 1.2) it is not. Never- 
theless, this limit law allows to expose the sample paths of £? mm (-). They are described 
below following [BF]. 

A generic trajectory of the process £? mm (-) is a continuous function and it is composed 
of decreasing, increasing and flat portions which follow each other in the following cyclic 
order: flat-increasing-flat-decreasing. A flat portion preceded by a decreasing portion and 
followed by an increasing one has length 1; such a portion is called valley. A flat portion 
preceded by an increasing portion and followed by a decreasing portion, has the length 9 
picked from the distribution P[6> < x\ = 2a a / 2 (l + x) _1 , x G [0, 1]. Such portion is called 
a plateau. An increasing portion is constructed from a realization of an increasing process 
C{x),x > 0, which is stopped at the height that has an exponential-1 law. The process £(■) 
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is a subordinator whose Levy measure \i (see [KaS] , Chapter 6) has the following form: 



= 2(^3)172 ' «e(0,l),d«c[0,l] 

A decreasing portion being taken with the sign "-" , is a probabilistic replica of an increasing 
portion. All the portions that compose I? mm (-) are independent which means that the 
random variables involved in the construction of these portions should be taken to be 
independent. 

We observe that what is lacking to have a complete description of £? mm (-) is the 
dependence between £? mm (0), the value of this process at x = 0, and the form of £? mm (-) 
as a whole. This happened because ( n preserves only the form of (M n fo)(-). We may 
present results in this respect in the future ([B]). 

Example 3.2. We shall use here our knowledge of the law of I? mm (-) in order to find 
approximately the distribution at time n of CA 184 with a particular initial distribution. 
Let rj n , n G N be a CA 184 whose initial state has the Bernoulli 1/2 distribution, that 

is 

i]o(i),i e Z, are i.i.d. with P [770(1) = l] = P [770(1) = 0] =1/2 (3.7) 

Let next g (-) denote the counting profile of 770 such that go(0) = 0. Then as it is well 
known, n-^goin-) -> B(-), where B(-) is as specified in Example 3.1 above. Thus, by 
Theorem 3.3, 

n- 1/2 g n (n-) - S min (-) (3.8) 

for an appropriate sequence of counting profiles. 

Our programme is now as follows: We shall base on (3.8) to assume that 

9n {.) = n^B^n- 1 -) in distribution (3.9) 

and shall combine this assumption with the law of B min (-) presented in Example 3.1 above, 
in order to derive the picture of the particle distribution in 77 n for large n. 

Although <7 n (-) is never a constant, n~ l / 2 g n (-) looks like a constant on an interval 
I C R in which either each odd site of 77 n is occupied by a particle while each even site is 
empty, or vice versa. Let us say that in this case we see a duce pattern of rj n in I. Similar 
considerations suggest us to call g n (-) increasing (resp., decreasing) on an interval I C R if 
77 n does not have a pair of contiguous sites of I occupied by particles (resp., holes) and at 
both endpoints one finds in r\ n a pair of contiguous sites occupied by holes (resp., particles) 
(the term "hole" used here, is a natural nickname for site free of particle) . Let us say that 
in this case we see a hole dominated (resp., particle dominated) pattern of r\ n in / 

Being equipped with the terminology of the above paragraph, we can now present the 
result of the programme we have layed down. It is the following picture of the distribution 
of 77 n : One will see the space R as divided in intervals in the way such that each one of 
them is either a particle dominated, or a hole dominated, or a duce pattern of 77^. The 
interval endpoints are distributed as homogeneous point process on R (whose parameters 
can be calculated explicitly from the law of £? mm (-)). The intervals follow each other in the 
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cyclic order: particle dominated - duce - hole dominated - duce - particle dominated. The 
mean extension of an interval is Cn + o(n), where C = C\ for both a particle dominated 
and a hole dominated regions, C = C2 for a duce region following a particle dominated 
region, and C = C3 for a duce region preceding a particle dominated region. The mean 
number of the pairs of contiguous sites occupied by particles (resp., free of particles) in a 
particle dominated (resp., hole dominated) region is C \/n + o(y / n). The numerical values 
of Ci, C 2 , C 3 , C may be calculated from the law of £? mm (-). 

Going back to the observation at the end of Example 3.1, we draw the reader's atten- 
tion to the fact that the information on the law of B min (-) that is absent did not affect the 
derivation of the law of r] n in the present example. The reason is in that the height of the 
counting profile g n (-) (i.e. <7n(0)) is indifferent to the distribution of the corresponding rj n . 

The picture presented above is not the exact law of r\ n because ^JnB Ulln {n~ l ■) might 
diverge from g n {-) by at most o(n 1 / 2 ) along the y-axis and by at most o(n) along the x-axis 
(in fact, they must diverge since the former does not attain its values in 1Z, and the latter 
does), and thus, (3.9) is an approximate relation. In fact, our argument does not control the 
particle amount with the precision higher than o(n 1 / 2 ), that is, the real amount of particles 
in a region / might differ from that of the picture presented, by at most |/| x o(n 1 / 2 ). Also, 
we were unable to give the particle distributions within each pattern, but instead, gave 
its characteristics in terms of its mean values. These are certainly, disadvantages of the 
presented approach when compared to a direct calculation of the law of rj n . An advantage 
is the simplicity: the presented approach requires only to know the law of W min (-) while a 
direct calculation requires a lot of work and produces usually cumbersome formulae. We 
note here that if one wishes to find i] n in the exact form then this could be done in the 
following way: translate CA 184 to BA (via T from Assertion 1.1), use then the tools 
developed in [BF] to find the law of BA at time n, and, finally, translate this law back into 
terms of CA 184. 

4. RATE OF CONVERGENCE TO THE EQUILIBRIUM 

In this section, we present an argument that estimates both the rate of convergence 
of CA 184 to its invariant state and the rate of the decay of particle density in BA (the 
former is the content of Theorem 4.1, and the latter is an auxiliary result established in 
the course of the proof of this theorem). The way this argument works is demonstrated 
here on a particular example which is the CA 184 starting from the Bernoulli 1/2 product 
measure. Remark 4.2 marks the place in this argument where main technical problems 
arise, when it is adapted to other initial distributions. 

Remark 4.1. Observe that that distribution at time n of the CA 184 studied here, has 
been approximately found in Example 3.2. From this approximation, one can conclude 
that the probability to see two fixed contiguous sites of Z at the same state in CA 184 at 
time n, decays at least as const x n -1 / 2 . A refinement of the arguments from Section 3 
would reveal that const x n -1 / 2 is in fact, the correct asymptotics for this decay. This 
is exactly the result of Theorem 4.1 presented below. We remark here that these two 
approaches have much in common, although apparently they are different. Indeed, if one 
possesses a tool to estimate (4.6) (which is the cornerstone of the proof of Theorem 4.1) 
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then the same tool may be employed to find the distribution of the process W min (-) which 
is the basic ingredient for the approximation argument exhibited in Section 3. 

Theorem 4.1. Consider CA 184 {r] n ,n G N} in which rjo is distributed by the Bernoulli 
1/2 product measure, i.e. {770(2), 7 G Z} are i.i.d. random variables with P [770(2) = l] = 
P [770(1) = 0] = 1/2. Let fi n denote the distribution of r\ n . Then (i n converges as n — > 00, 
to 1/2 (<5 + S e ) at the rate n -1 / 2 . In exact terms this means that for each k G N there 
exist < c' k < Ck < 00 such that 

c' k n- 1/2 < sup \fi n (U)-l/2(5 + 5 e )(U)\<c k n- 1 / 2 (4.2) 
c/e{o,i}{°. 1 . ■■■■'=} 

Proof. Fix k G N and define if := {0, 1, . . . , fc}. Define 

-E 1 := {77 G {0,1}^ : 3 i G {0, 1, . . . , k - 1} for which 77(1) = 77(2 + 1)} 

Let a, (3 G {0, 1} K be the restrictions of respectively, o and e (the configuration e and o 
have been defined in Theorem 2.2) to K : a(0) := 0, (3(0) := 1, a(i + 1) := 1 — 0(7), 
+ 1) := 1 - P(i), i = 0, 1, . . . , k - 1. Then 

i?U{a,/3} = {0,l} K (4.3) 

We now observe that since /io is translation invariant then /j n must be such as well, and 
this fact leads to 

A»n(a) = /«n(/3) (4.4) 

Since l/2(5 + 5 e )(E) = and l/2(<5 + <5 e )(a) = 1/2(S + S e )((3) = 1/2 then (4.3) and (4.4) 
imply that the supremum in (4.2) is fi n (E). In order to estimate this measure we shall use 
BA {Cn,n G N} such that 

Co = Ti] where rj ~ Bernoulli 1/2 product measure and T is defined by (1.3) (4.5) 

^From the definition of Co above and Assertion 1.1, we have that 

Hn(E) = P [CnW + for some % G {1, 2, ... , fc}] 

Let L = L(n) (R = R(n)) denote the event that (in the considered BA) a particle 
that originated from {— n, —n + 1, . . . , — n + k} (resp., {n, n + 1, . . . , n + k}) is in the region 
K' := {1,2,..., k} at time n. Obviously, 

P[L] < fi n (E) =P[Ln.R c ] +P[L c n.R] +P[Ln.R] < P[L] +P[R] = 2P[L] 

whereas to derive the last inequality we used that changing Co to —Co does not affect the 
distribution of Co defined in (4.5). As for P[£], one easily gets the following quite rough 
bounds 

cP < P \L\ < ckP, for an appropriate < c < 00 
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where P denotes the probability to have a positive particle at site at time and for this 
particle not to have survived till time n (observe, this is the same value, if is substituted 
by any site of Z). 

To estimate the value of P we we will use the following assertions. Their proofs are 
straightforward. 

Assertion 4.1. A particle in the BA will be annihilated till time n if and only if the 
distance between its initial position and that of its annihilating companion is < 2n — 1. 

Assertion 4.2. Let Co be a configuration from {—1,0, 1} Z . Assume ( (i) = 1 for some 
i. Let /(•) be the integrated profile corresponding to Co and such that f(i) = 0. Then 
min{/c > i + 2 : f(k) = 0} — 1 is the initial position of the annihilating companion in BA 
{Cn, n £ N} of the particle that started from i. 

Let /(•) be the integrated process related to Co whose distribution is (4.5) conditioned 
to having a positive particle at 0. Then due to Assertions 4.1 and and 4.2, we conclude 
that 

P = P [/(0) = 0, /(l) > 0, . . . , f(2n - 1) > | /(l) > 0] (4.6) 

There are several ways to get the decay rate of (4.6). The promptest one goes along the 
following lines: The random variables {Co(*)}iez are not independent, but fall under the 
conditions of Chapter 4 of [Bi], which provide that cn -1 / 2 f(n-) — > B(-) for some constant 
c. From this one concludes that (4.6) ~ const x n -1 / 2 . This completes the proof. A 

Remark 4.2. The applicability of the presented proof to another distribution of Co (and 
thus, also of 770) depends basically on ones ability to calculate the asymptotics of (4.6) for 
this distribution. 
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